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Abstract 

' We start with a simple introduction into the renormahzation group (RG) in 

■ quantum field theory and give an overview of the renormalization group method. 
, The third section is devoted to essential topics of the renorm-group use in the 

0\ ■ QFT. Here, some fresh results are included. 

■ Then we turn to the remarkable proliferation of the RG ideas into various fields 
, of physics. The last section summarizes an impressive recent progress of the "QFT 

renormalization group" application in mathematical physics. 
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1 Renormalization group primer 



1.1 Mathematical preliminaries 

1.1.1 Renorm-group folklore 

Let us start with some simple statements which can be supposed to be widely known by 
particle theorists. In the quantum field theory (QFT) the renormalization group (RG) is 
usually associated with a possibility of presenting any physical quantity, F{Q'^,g), cal- 
culated under a definite renormalization prescription in the form F [Q"^ / ^'^ , g ^) (for 
simplicity - in a massless case) with the renormalized coupling constant (7^ definition 
attached to some renormalization point (or reference momentum scale) Q = fi. Differen- 
tial RG equation is usually said to be driven from the condition that F does not depend 
on the choice of fj, , 

dF 

The coupling constant dependence on /x is described by a specific function g{Q'^) 
known as an effective coupling (sometimes - effective coupling constant) g^ = ^(yU^) . 
Eq. (0) can be written down in the form of a partial linear differential equation (DE) 



dx ^^^^ dg 



F{x,g) = 



(2) 



where x = Q'^/fi^, g stands for g^^ and P{g), the group generator, usually referred to as 
beta function, is defined by 



Pig,) 



dg{z) 
dz 



at z = fx . 



The effective coupling g should be considered as a function of two arguments: x 
Q'^/fx^ and g^ with the boundary condition g{l,g) = g. Besides 



dx dg 



9ix,g) = 



(3) 



it satisfies the nonlinear DE 



X- 



d g{x,g) 
dx 



Pi9ix,g)) 



(4) 



which is nothing else but a characteristic equation for (D). To employ this formalism, 
one has to give l3{g). Usually, for this one uses renormalized perturbation theory. 

The foregoing can be considered as a "RG folklore". For brevity, we gave it in the 
simplest massless version, which corresponds to the UV case, for the QFT model with 
one coupling constant. 
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1.1.2 Group Functional Equation 

Less popular are the RG Functional Equations (FEs). The FE for the 'g in the UV case 
has the form 

9{x,g) = g(^j, g{t^9)) ■ (5) 

This equation, which follows (see, e.g., the Chapter Renormalization Group in Ref.[jl|) 
from finite Dyson renormalization transformations, represents a basement of the differen- 
tial RG formulation. Popular DE (^) can be directly obtained from it by differentiating 
over X and then putting t = x. On the other hand, by differentiating (^) with respect 
to t at if: = 1 we get partial DE (|^). 

The FE (^) as well as similar FEs for propagators and vertex functions (see, below, 
eq.(|T^) must be considered as the most compact and general formulation of the RG 
symmetry in QFT. 

However, in reality, group FEs, like (^) (and DEs (^ and (^ as well) do not contain 
any physics at all being just the reflection of the group composition law! Here, we mean 
the continuous group (that is, the Lie group of transformations) of operations changing 
the reference point /i involved into the coupling constant (yf^ definition. Namely, we can 
regard the change of a reference coupling g^ as an operation of the group element 
Tt 

Ttg^ = g^^t = g{t, g,,) 

with a real continuous positive numerical parameter t (= /i^/ /2^) . 

If we set X = rt, then the l.h.s. of (|^) can be achieved from g by operation Trt, while 
the r.h.s. may be identified as Tj-Ttg. The content of eq.(||) is just the group composition 
law, 

Trt = TrTt . 

Thus, the essence of the basic RG functional equation (|) is the necessary condition 
for transformations Tt to form a group. 

At the same time, it demonstrates that function g is invariant with respect to simul- 
taneos transformation 

Rt : {x' = x/t, g' = g{t,g)} . (6) 
Invariance condition for an observable now can be written down as 

F{x,g) = F(^j, git,g)^ . 

Usually, of interest are also functions (f){x,g) (like, e.g., propagator amplitudes in 
QFT) transforming as a linear representation of RG 

(f){x,g) Rt(f) = (f){x,g') = z{t,g)(f){x,g) . (7) 

Note also that the group FE for an observable, like matrix element, is of the form 

Mi{x},y;g)=M(^S^jYj^9{'t,y;g)J; {x} =Xi,X2,...,Xk (8) 

which reflects an existence of several Q^-type arguments and implements its indepen- 
dence of renormalization details corresponding to Eq.([l|). 
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1.1.3 Abstract Formulation 



To make this point clearer, let us show that FE can formally be obtained directly 
from the group composition law. Generally, the mathematical formulation of the RG 
transformation can be presented as a functional realization of the mentioned Lie group. 

Consider the transformation T(/) of a certain abstract set A4 of elements Mj into 
itself depending on a continuous real parameter / (—00 < / < 00) such that for each 
element M we have 

T(/)M = M' {M,M' C M) . 

Suppose that M. can be projected onto the real axis, i.e., for every Mj there corresponds 
a real number g^. Then, this transformation can be written in the analytic form 

T{l)g = g' = G{l,g) , 

G being a continuous function of two arguments satisfying the normalization condition 
G(0, g) = g which corresponds to the identity transformation T(0) = E. 
Transformations T{1) form a group if they satisfy the composition law 

T(0 X T(A) = T(/ + A) (9) 

to which there corresponds the functional equation for G : 

G{l,GiX,g)} = Gil + X,g) . (10) 

As it follows from the bases of the Lie group theory, it is sufficient to deal with the 
infinitesimal transformation at A -C 1 , i.e., with the DE 

^^ = /3{G((.9)}. (11) 

Here the group generator is defined as 

dG{e,g) 
I3[g) = — — at e = 0. 

Performing a logarithmic change of variables 

l = \nx, X = lnt, G{l,g) = g{x,g) , T{lnt)=Tt (12) 

we obtain (|^) and (^ instead of (0) and (|1T]). 

-'^This condition is not essential and can be modified — see, below, eqs. (|6|) and (p7|). 
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1.2 Definition of the Renorm-Group 

1.2.1 The RG transformation 

Generally, the RG can be defined as a continuous one-parameter group of specific trans- 
formations of a partial solution (or the solution characteristic) of a problem, a solution 
that is fixed by a boundary condition. The RG transformation involves boundary condi- 
tion parameters and corresponds to some change in the way of imposing this condition. 

For illustration, imagine an one-argument solution characteristic f{x) that has to 
be specified by the boundary condition /(xq) = /o- Formally, represent the given char- 
acteristic of a partial solution as a function of boundary parameters as well: f{x) = 
f{x,xo,fo)- (This step can be considered as an embedding operation). The RG trans- 
formation then corresponds to a changeover of the way of parameterization, say from 
{^0, fo} to {xi, fi} for the same solution. In other words, the x argument value, at which 
the boundary condition is given, does not need to be Xq, but we may choose another 
point Xi. Our solution / can be written in a form of a two-argument function F{x/xo, fo) 
with the property F(l,7) = 7. The equality F{x/xo, fo) = F{x/xi, fi) reflects the fact 
that under such a change of a boundary condition the form of function F itself is not 
modified (as, e.g., in the case of F(x, 7) = $(lnx + 7)). Noting that fi = F{xi/xo, fo) , 
we obtain 

^(e, fo) = Fi^/t, F(t, fo)) ; e = x/xo , t = xjxo . 

The group transformation here is { ^ ^ ^/t, fo ^ fi = F{t, fo) } . 

The renorm-group transformation for a given solution of some physical problem in the 
simplest case can be defined as 

a simultaneous one-parameter transformation of two variables, say x and g, by 

Rt : {x^x' = x/t, g ^ g' = g{x,g) } , (6) 

the first being a scaling of a coordinate x and the second — a more complicated functional 
transformation of the solution characteristics. Eq.(^ for the transformation function 'g 
provides the group property T^^ = T^-Tt of the transformation (6). Performing the loga- 
rithmic change of variables and an appropriate redefinition of a transformation function 
(|13), we obtain eqs.(|ig), (|l|) and 

m) ■■ {q^q' = q-l, g ^ g' = 0(1, g) } , (13) 

instead of (|), (^) and (6). One can refer to these equations as the multiplicative version 

(and previous equations in abstract formulation as to the additive one). They are just the 
RG equations and transformation for a massless QFT model with one coupling constant. 
In that case x = is the ratio of a 4-momentum Q squared to a "normalization" 

momentum squared and g, the coupling constant. 

Several generahzations of @ and ([13D will be considered below. 
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1.2.2 Simple generalizations 



"Massive" Case 

For example in QFT, if we do not neglect the mass m of a particle, we have to insert 
an additional dimensionless argument into the invariant coupling g which now has to be 
considered as a function of three variables: x = Q'^ j jj?, y = w? / jj?, and g. The presence 
of a new "mass" argument y modifies the group transformation 

Rt : {a;' = | , y' = | , g' = -g{t.y:9) ] (14) 

and the functional equation 

9{x,y;g) = g , - ] g{t,y]g)j . (15) 

Here, it is important that the new parameter y (which in physical nature must be close 
to the X variable as it scales similarly) enters also into the transformation law of g . 

If the considered QFT model, like QCD, contains several masses there will be several 
mass arguments 

y ^ {y} = yuy2,---yn ■ 



Multi-coupling case 

A more complicated generalization corresponds to transition to the case with several 
coupling constants: g {g} = gi, . . . gk . Here, one has to introduce the "family" of 
effective couplings 

9 ^ {9} , 9i = 9ii^,y;{9}) ; z = i,2, ...fc, 

satisfying the system of coupled functional equations 

f X y \ 

9i{x, y] {9}) = 9i{-, - ; { 9{t, y, {9}) } j ■ (16) 

In the abstract formulation this system is a generalization of (^ and ([15|) for the case 
when every element Mi of Ai can be described by k numerical parameters, i.e., by a 
point {g} in the fc-dimensional real parameter space. The RG transformation now is 

Rt ■■ {^-^7' {9}^{9m] , 9^{t)=Ut,y■A9}) ■ (17) 

1.3 Early history and the RG method 

1.3.1 Renormalization and renormalization invariance 

As it is known, the regular formalism for eliminating the UV divergences in QFT was 
developed on the basis of covariant perturbation theory for the scattering S'-matrix in 
the late 40s. This breakthrough is connected with the names of Tomonaga, Feynman, 
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Schwinger and some others. In particular, Dyson and Abdus Salam carried out the 
general analysis of the structure of divergences in arbitrarily high orders of perturbation 
theory. Nevertheless, a number of subtle questions concerning overlapping divergences 
remained unclear. 

An important contribution in this direction based on a thorough analysis of the 
mathematical nature of UV divergences was made by Bogoliubov. This was achieved on 
the basis of a branch of mathematics which was new at that time, namely, the Sobolev- 
Schwartz theory of distributions. The point is, that propagators in local QFT are 
distributions (similar to the Dirac delta-function) and propagator products appearing 
in the coefficients of the S'-matrix expansion require a supplementary definition in the 
case when their arguments coincide and lie on the light cone. 

In the mid 50s on the basis of this approach Bogoliubov and his disciples developed a 
technique of supplementing the definition of products of singular Stiickelberg-Feynman 
propagators [|| and proved a theorem 0] on the finiteness and uniqueness (for renor- 
malizable theories) of the S'-matrix elements in any order of perturbation theory. The 
prescriptive part of this theorem, the Bogoliubov R-operation (see, e.g., chapter "Removal 
of divergencies from the S'-matrix" in the monograph |I|), still remains a practical means 
of obtaining finite and unique results in the higher order perturbation calculation. 

The Bogoliubov algorithm works, essentially, as follows: 

- To remove the UV divergences of a one-loop diagram, instead of introducing some 
regularization, e.g., the momentum cutoff, and handling (quasi) infinite counterterms, it 
suffices to complete the definition of a divergent Feynman integral by subtracting from 
it a certain polynomial in the external momenta which in the simplest case is reduced 
to the first few terms of the Taylor series. 

- For multi-loop diagram (including one with overlapping divergences) one should 
first subtract all divergent subdiagrams and finish with subtracting the diagram as a 
whole in a hierarchical order regulated by the i?-operator. 

An attractive feature of this approach is that it is free from any auxiliary nonphysical 
attributes such as bare masses, bare coupling constants, and regularization parameters 
which are not involved in the computation within the Bogoliubov's algorithm. The latter 
can be regarded as renormalization without regularization and counterterms. 

The uniqueness of computational results for the observable S'-matrix elements is 
ensured by special conditions imposed on them. These conditions contain some degree 
of freedom (related to different renormalization schemes and momentum scales) that 
can be used to establish finite relations between the Lagrangian parameters (masses, 
coupling constants) and corresponding renormalized quantities. The fact that physical 
predictions are independent of arbitrariness in the renormalization conditions, that is, 
they are renorm-invariant, constitutes the conceptual foundation of the renormalization 
group. 
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1.3.2 The discovery of the renormalization group 

In the 1952-1953 Stiickelberg and Peterman discoveredQa group of infinitesimal trans- 
formations related to finite arbitrariness arising in the S'-matrix elements upon elimina- 
tion of the UV divergences. These authors introduced normalization group generated by 
Lie operators connected with renormalization of the coupling constant e. 

In the following year, on the basis of (infinite) Dyson's renormalization transfor- 
mations formulated in the regularized form, Gell-Mann and Low derived functional 
equations for the QED propagators in the UV limit. The appendix to this article con- 
tains the general solution (obtained by T.D. Lee) of this functional equation for the 
renormalized transverse photon propagator amplitude (i(x, e^), written in two equiva- 
lent forms: 



e^d(x,e^) = F {xP-^ {e^)) and Inx = J —fr ^ ^(e 



^Pdl) d\nx 



x=l 



A qualitative analysis of the behaviour of the quantum electromagnetic interaction at 
small distances was carried out with the aid of (p^). Two possibilities, namely, infinite 
and finite charge renormalizations were pointed out. 

However, paper paid no attention to the group character of the analysis and the 
results obtained there. The authors missed a chance to establish a connection between 
their results and the standard perturbation theory and did not discuss the possibility 
that a ghost pole solution might exist. 

The final step was taken by Bogoliubov and the present author bs-55a,bs-55b,sh-55 
— see also the survey [|10| published in English in 1956. Using the group properties of 
finite Dyson transformations for the coupling constant, fields and Green functions, these 
authors derived functional group equations for the propagators and vertices in QED in 
the general case (that is, with the electron mass taken into account). For example, the 
equation for the transverse amplitude of the photon propagator and electron propagator 
amplitude were obtained in the form 

d{x, y; e^) = d{t, y; e^) d (^^ , |; e^d{t, y; e^) j , 

s{x, y; e^) = s{t, y; e^) s (^^, |; e'^d{t, y- e^)^ (19) 

in which the dependence on the mass variable y = m? / was present. 

As can be seen, the product e^d of electron charge squared and photon propagator 
amplitude enters in both the EEs. This product is invariant with respect to Dyson's 
transformation. We called this function - invariant charge and introduced the term 
renormalization group. 



For a more detailed exposition of the RG early history see our review 
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In the modern notation, the first equation is that for the invariant charge (now widely 
known as an effective or running coupling) d = ad{x, y;a = e^): 

a{x,y-a) = a{-,j] a{t,y;a)] . (20) 

Let us emphasize that, unlike the approach Ref. [|], in the latter case there is no 
relation with UV divergences and simplification due to the massless nature of the UV 
asymptotics. Here, the homogeneity of the transfer momentum scale is violated explic- 
itly by mass m. Nevertheless, the symmetry (even though a bit more complex one) 
underlying the RG, as before, can be stated as an exact symmetry of the solutions of the 
QFT problem. This is what we mean when using the term Bogoliubov's renormalization 
group or renorm-group for short. 

The differential group equations for a and for the electron propagator: 
da{x,y]a) (y _ \ ds{x,y;a) (y _ \ 

d\^x = ^ U' ^' "V ' = ^ U' ^' "V ^' ' ^''^ 

with 

da{^,y;a) ds{^,y;a) 

P{y^(^) = . 7(2/, a) = at^ = l (22) 

were first derived in by differentiating the FEs. In this way, explicit realization of 
the group DEs mentioned in the paper was obtained. These results established a 
conceptual link with the Stiickelberg — Peterman and Gell-Mann — Low results. 

1.3.3 Creation of the RG method 

Another important achievement of paper consisted in formulating a simple algorithm 
for improving an approximate perturbative solution by combining it with the Lie equa- 
tions (for detail, see below. Section 2). 

In our adjacent publication |p this algorithm was effectively used to analyse the UV 
and infrared (IR) behaviour in QED. The one-loop and two-loop UV asymptotics 



-(1) I \ — -W t n \ ^ 

x'r^[x- a) = a)j^(x, 0, a) = ^ _ ^^^^ 



(23) 



3n 



,,(2) . . a 



^ l-flnx + |fln(l-tlnx) ^''^ 
of the photon propagator as well as the IR behavior 

s{x, y; a) ^ {x/y - l)-W2- = (p^m^ - ly^"/^- (25) 

of the electron propagator in the transverse gauge were obtained. At that time, these 
expressions were already known only at the one-loop level. It should be noted that in 
the mid 50s the problem of the UV behaviour in local QFT was quite urgent. At that 
time, substantial progress in the analysis of QED at small distances was made by Landau 
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and his collaborators |TT[|. However, Landau's approach did not provide a prescription 



for constructing subsequent approximations. 

The simple technique for obtaining higher approximations was found only within the 
new renorm-group method. The one-loop UV asymptotics of QED propagators obtained 
in our paper eqs. (|23| ) and (pSD , agreed precisely with the results of Landau's group. 



Within the RG approach these results can be obtained in just a few lines of argu- 
mentation. To this end, e.g., the massless one- loop perturbation approximation should 
be substituted into the r.h.s. of the first equation in ( P^ to compute the generator 
/3(0, a) = ipi^a) = a^/Svr followed by elementary integration of the first of eqs.(^). 

Moreover, starting from the next order perturbation expression apj;,f^{x, ; a) con- 
taining the In x term, we arrived at the second renorm-group approximation (p^ ) 
performing summation of the a^(aln)"' terms. This two-loop solution for the invariant 
coupling first obtained in contains the nontrivial log-of-log dependence which is now 
widely known of the "next-to-leadind logs" approximation for the running coupling in 
quantum chromodynamics (QCD) — see, below, eq . (^Tf) . 



Comparing solution (^) with (p3D, one can conclude that the two- loop correction is 
extremely essential just in the vicinity of the ghost pole singularity at xi = exp {3n/a). 
This demonstrates that the RG method is a regular procedure, within which it is quite 
easy to estimate the range of applicability of the results. 

Quite soon, this approach was formulated for the case of QFT with two coupling 
constants, say, g and h, namely, for a model of pion-nucleon interactions with the pions 
self-interaction. To the system of functional equations for two invariant couplings 

f{x,y; g^,h)j=f(^^, |, f{t,y;g^,h), h(t,y; g'^.h^^ , 
h (x, y; g^, =h (^j, |, f (t, y; g^, /i) , h (t, y; g^, h 

there corresponds a coupled system of nonlinear DEs - see, below, eqs.(^). It was 
analysed [0 in the one-loop appriximation to carry out the UV analysis of the renor- 



malizable model of the pion-nucleon interaction. 

In a more general case of arbitrary covariant gauge the RG analysis in QED was 
carried out in Here, the point was that the charge renormalization is connected 



only with the transverse part of the photon propagator. Therefore, under nontransverse 
covariant (e.g., Feynman) gauge the Dyson transformation has a more complex form. 
This issue has been resolved by considering the gauge parameter as another coupling 
constant. 

In Refs.[0, H, |13| and [l^ the RG was thus directly connected with practical 



computations of the UV and IR asymptotics. Since then this technique, known as the 
renormalization group method (RGM) and being summarized in the first edition of 
monograph M, has become the sole means of asymptotic analysis in local QFT. 



1.4 RG in QED 
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1.4.1 Effective Electron Charge 

An essential feature of quantum theory is the presence of virtual states and transitions. 
In QED, e.g., the process of virtual dissociation of a photon into an electron-positron 
pair and vice versa 7 ^ e"*" + e~ can take place. The sequence of two such virtuale 
transitions represents the simplest contribution to the effect of vacuum polarization. 

The vacuum polarization processes lead to several specific phenomena and partic- 
ularly to the notion of effective electron charge. To explain this, let us start with a 
classical analogy. 

Take a polarizable medium consisting of molecules that can be imagined as electric 
dipoles. Insert into it an external electric charge Q. Due to the attraction of opposite 
charges, the dipoles change their position so that the charge Q turns out to be partially 
screened. As a result, at a distance r from Q the electric potential will be smaller than 
the vacuum Coulomb law Q/r and can be presented in the form Q{r)/r where, generally, 
Q{r) < Q. The introduced quantity Q{r) is known as an effective charge. As r decreases, 
Q{r) increases and as r — 0, Q{r) tends to Q. 

In QFT the vacuum, i.e., the interparticle space itself stands for the "polarizable 
medium". Quantum-field vacuum is not physically empty. It is filled with vacuum 
fluctuations, i.e., with virtual particles. These "zero fluctuations" are a well-known 
effect of a ground state in quantum world. In QED, zero oscillations consist mainly of 
short-lived virtual {e~^,e~) pairs which play the role of tiny electric dipoles. 

Consider the process of measuring the electron charge with the help of some external 
electromagnetic field. In the quantum case the probing photon can virtually dissoci- 
ate into the (e+,e~) pair. This pair can be treated as a virtual dipole that produces 
partial screening of the measured charge. The simplest process involves two elementary 
electromagnetic interactions, its contribution to an effective charge being proportional 
to the small number = a ~ 1/137; and this contribution depends on the distance 
r ! In the region of r values much smaller than the Compton length of the electron 
= h/mc ~ 3,9.10~^^cm it depends on r logatithmically 

e e(r) = e(l-— In— + (26) 
I 3vr J 



as was first discussed by Dirac in the middle of the 30s. The e(r) value decreases as 



r grows. So, qualitatively, the QED effective charge behavior corresponds to a classical 
picture of screening. 

This dependence can be presented by a set of curves e(r). Each curve represents a 
possible behavior of the effective charge e(r) as obtained from the theory and considered 
without any reference to experiment [a = being unspecified numerically). 

The point is that in the classical analogue the value of an external charge Q inserted 
into the polarizable medium is known from the very beginning. In quantum physics it 
is not the case, and a charge value can be measured at not very small distances. The 
result of measurement generally has to be specified by two quantities: the "distance of 
measurement" and the measured charge value Cj. Hence, to make the choice from 
the mentioned set of curves, one has to fix the point on the plane with the coordinates 
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r = Tj, e(r) = Ci. Thus, for the chosen "physical" curve e(rj) = Cj. Note, that the usual 
definition of the electron charge by a classical macroscopic (like Millikan) experiment 
corresponds here to very large distances r > Vg, i.e., e = e(r = re) = l/vl37. 

As it is well known, in relativistic microphysics one usually uses the momentum rather 
than coordinate representation. Correspondingly, instead of e(r) one deals with the 
quantity the Fourier transform of e(r) squared. It is a monotonically increasing 

function of its argument Q^, the 4- momentum transfer squared. Here and below, the bar 
denotes a function (distinct from a, a^, ctj - its numerical values at some given value of 
the argument). The correspondence condition with the classical electrodynamics now 
takes the form a(0) = 1/137, as in our scale the external long-range field corresponds to 
a photon with vanishing 4-momentum. However, as before, to fix one of possible curves 
on the plane {Q, a) one has to give a point Q = ^fQ^ = /i, a = and hence, for the 



a 



selected curve Q;(/i^) 

The parameter /x sometimes is referred to as a scale 'parameter. As is clear, it is just 
the momentum magnitude for a photon used for the charge measurement. The effective 
coupling function a((5^) describes the dependence of the electron charge value on the 
measurement conditions. In our days the logarithmic corrections to the Millikan value 
become essential and are measured at big accelerators. 

The parameter [i has no analogue in the QED Lagrangian that reproduces the clas- 
sical electro dynamical one. The phenomenon of its arising in QFT was exaggerated by 
the term "dimensional transmutation" . As it was shown, its appearance is very natural 
and is connected with the measurement procedure. 

This is a good place to recall the ideas by Niels Bohr formulated in the middle of 
30s [0 and related to the complementarity principle. The point is that to specify a 
quantum system, it is necessary to fix its "macroscopic surrounding", i.e., to give the 
properties of macroscopic devices used in the measurement process. Just these devices 
are described by additional parameters, like /i. However, this is not the end of the Bohr 
(i.e. - scale) parameter story. As can be shown, in the QFT this parameter existence 
leads to a new symmetry lying in the foundation of the renormalization group. 



1.4.2 RG transformations 

To do this, consider again the function a{Q'^) having in mind that the physical solution 
has been chosen by the condition a(Q^ = /i^) = a^. Assume also for simplicity that 
we deal with a massless QED, more precisely, with the approximation \ Q \^ m. This 
corresponds to the GeV-energy region or to distances r ^r^,. Here the effective charge 
function a can be represented as a function of two dimensionless arguments Q"^ / /i^ and 
a^, i.e., a(Q^) = a{Q'^/^^,a^) . 

Now, take into account that the couple of parameters used to identify the 

physical solution may, generally, correspond to any scale /x. Take two scales "1" and "2" 
with coordinates /xi, ai and /i2, ^2, respectively. It is evident that a can be parameterized 
by any pair /Xj, a^; z = 1, 2, . . . so that for arbitrary values the identity 

«(QVAii5"i) = a(QV/"25"2) 
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should hold. 

At the same time the second argument in the r.h.s., a2, which by definition is equal 
to a at = /U2, can be expressed in terms of a parameterized with the help of point 
"1" coordinates, i.e., 02 = = H2) = (^{lAl lAi ^i) ■ Combining the last two relations 
and introducing a notation j ^ = x, a\ = a, = t , we arrive at the FE 



identical with eq.(^. 

Note, that the corresponding continuous one-parameter transformation is just the 
change (" 1" " 2" ) of the parameterization point 



Thus, we have shown that in the renormalized QED there exists invariance with 
respect to continuous transformations of the group type which involve two quantities 
and contain a functional dependence. This precisely corresponds to definition (^. 

As can be shown, in QED the effective coupling a is equal to a product of a and 
dimensionless function d{x, a) - the transverse photon propagator amplitude with due 
regard for vacuum polarization effects. Generally, in QFT models with one coupling 
constant the invariant coupling g{x, g) can be expressed as a product of g, corresponding 
vertex function and the square root of propagator amphtudes of the fields participating 
in the interaction. Usually, this can be done on the basis of Dyson finite renormalization 
transformations. 

Thus, the RG invariance is nothing else but the invariance of a solution with respect 
to the way of its parameterization. For instance, in real QED, instead of using the 
"MiUikan's value" a(0) = 1/137 one may take the "CERN value" a(M|) ~ 1/128, 9. 

2 Renormalization group method 
2.1 Basic idea 

Approximate solution of the physical problems with RG symmetry usually does not 
obey this symmetry which is lost in the course of approximation. This is essential when 
the solution under consideration posseses a singularity as far as the singularity structure 
commonly is destroyed by an approximation. 

In QFT, e.g., the usual way of calculation is based on the perturbation method, i.e., 
on power expansion in ^f. It is not difficult to see that finite sums of this expansion 
do not satisfy the functional group equations. As the simplest illustration, consider the 
effective coupling g in the UV region where the one-loop contribution has a logarithmic 
form 




(27) 



Rt '■ ^ = VijU^i, «! ^ ^2 = C({t, Ctl)} 





(29) 
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with /3, a numerical coefficient. By substituting this expression into FE after simple 
manipulation one has 

Discr[g^^] = gi^-l{x,g) - (^^t (a^A, ^7pt(^> fl')) = 

= [g + g'^plnx] - [g + g'^plnx + 2g^(3'^lntln{x/t)] ^ 

— error in the g^ order. This discrepancy can be liquidated by addition of the particular 
next order term to the r.h.s. of 



g^^l = g + g^f3\nx + g^(3^ In^ x. 

This "improved" expression would yield the discrepancy of the g"^ order which in its 
turn can be abolished by adding the g^ In^ term to ( ppj ) and so on. 

Thus, we see, on the one hand, that the finite polynomials cannot satisfy the condition 
of renormalization invariance. On the other hand, we can conclude that the functional 
RG equation represents a tool for iterative reconstruction of renorm-invariant expression 
that has the form of infinite series. 

This example illustrates a rather general situation. As a rule, approximate solutions 
do not satisfy a group symmetry. Here, this happens in the UV limit as In x ^ oo where 
the observed discrepancy becomes important. 

Another illustration is provided by the one-dimensional transfer problem (for detail, 
see our reviews in Refs. |T^). Take a half-space (/ > 0) filled with a homogenious media. 



Let some given amount of particles (or radiation) be falling on the surface (at / = 0) 
from the empty half-space. The particle density n{l,\) is a function of coordinate and 
of particle velocity v. It satisfies an integro-differential kinetic Boltzmann equation. 
In some cases one can neglect by the energy dependence of cross-sections. Here, the 
solution can be treated as a function of coordinate and direction of particle velocity 
= v/t>. In this, "one- velocity" , case the simple symmetry of the RG type was found 
not for density, but for n integrated over directions in forward hemisphere 



G{1) = f n{i,n)dn 



The function G relates to amount of all particles moving inwards the media. A par- 
tial solution of this problem will depend on the boundary condition at / = 0. Corre- 
spondingly, the solution characteristic G will be the function of two arguments G{l,g) 
- coordinate / (distance from the boundary) and total amount of ingoing particles 
g = G{1 = 0) = G{0,g). Just this function G{l,g) satisfies the group FE ( pUj ) 
in the additive form. 

It is rather simple to get an approximate behavior of this density G{l,g), at small / 
G{l,g) = g + lG\0,g), / « 1, (30) 

which, being considered for large / values, also does not obey the mentioned symmetry. 
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On this basis one can set the task of "renormahzation-invariant improvement" of 
perturbative results. The key idea is to combine an approximate solution with the 
group equations. The simplest and most convenient way for this "marriage" is the use 
of Lie equations, i.e., group differential equations. The renormalization group method 
(RGM) as it was first formulated in Refs.0, ||, |l^ is essentially based on these group 
equations. 



2.2 Differential formulation 



The differential equations can be obtained from the functional ones in two different ways. 
Differentiating eq.(|T5|) by x and putting then t = x , one obtains (compare with (H)): 



X- 



dg {x,y,g) 
dx 



dg{t,y;g) 



dt 



(22a) 



t=i 



The nonlinear equation (22a) can be considered as "massive" generalization of eq.(p. 
On the other hand, one can differentiate eq.(|T3[) with respect to t at the point t = 1, 
which yields 

'''^ + y'^~f^^y^3)^9{x,y;g) = 0, (31) 

a linear partial differential equation (PDE). 

Analogous operations applied to the second of eqs.(|l9D lead to: 



ds{x,y;g) 
d Inx 



7 



and 



where 



y ( 
-,9{x,y; 

X 



d 



s{x,y;g) 



d d 



dx 



i{y,9) 



ds{t,y]g 



dt 



(226) 
(32) 
(33) 



i=l 



is the so-called anomalous dimension of s. For a group invariant, like, e.g., matrix 
element M satisfying FE (|^) this dimension is equal to zero. The corresponding PDE 
looks like 



(34) 



Equations (|3TD , (pa), and (|3^) express the independence on the t parameter of the 
r.h.s. of the related functional group equations, i.e., a mutual compensation of t depen- 
dences via three (or more) arguments. This DEs can be called compensational equations 
to distinguish them from nonlinear eqs.(plD which can be referred to as evolutional group 
equations. 

Stress that compensational as well as evolutional DEs taken together with normaliza- 
tion (i.e. boundary) conditions like g{l,g) = g , s{l,g) = 1 are equivalent to functional 
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equations and to each other. At the same time, evolutional Lie equations turn out to be 
more convenient for practical construction of the solution, generators /3, 7 being given. 
Let us comment also that the UV limit of compensational DEs like, e.g., 

[^^-(3{9)§-g-l{9)^s{x,g) = (35) 

coincides with the UV limit of specific nonclosed equation 

- f3{g)-^ - 7(^?)| s{x, g) = AS {C - S) 

obtained in the early 70s by Callan and Symansik. The r.h.s. of this equation contains 
the result of mass counter-term insertion into all internal lines of all diagrams for the 
function s under consideration. For this reason, in current literature compensational 
equations are often related to as the Callan-Symansik equations. However, these equa- 
tions just in the form (0) and (|32|) were first obtained by Lev Ovsyannikov in 1956 



while solving ||T^ functional RG equations. Therefore, we consider it justifiable to relate 
compensational DEs to the Ovsyannikov's rather than to some other names. 

It is not difficult to formulate group DEs for a multi-coupling case by proper differ- 
entiation of EEs ([161) . For instance, the system of evolutional DEs looks like 

d9i{x,y,{g}) _ a fy 



dx 



A(^,te(x,y;M)}) . (36) 



2.3 General solution 

General solution of the group FEs was obtained in the paper ||18| by applying the theory 
of PDE to the compensational eqs. (|3T|) and (^). Details of the derivation can be found 
in the Section 48.3 of the third edition of the monograph |jl[|. The results obtained can 
be formulated as follows: 

To every solution of DE (0), there corresponds some function of two arguments 
F{y^ g), reversible with respect to its second argument and connected to g by the relation 

F{y,g) = Fi^^, g{x,y;g)^ . (37) 

The explicit form of g can be obtained now by reversing the r.h.s: 

9{x, y, g) = F^^ 

To determine F it is sufficient to specify the generator (3{y,g). 

Note also that to get from the Ovsyannikov result ( ^7\) the solution in the UV limit, 
i.e., in a massless case at y = 0, one has to assume for F a specific limiting form 

Fiy,g) = yexp[f{g)] or =\ny + f{g) as y^O. 



-, F{y,g) 

.X 
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Then 

f{-g{x,g)}-f{g) = \nx; g = f-^{\nx + f{g)} . (38) 

Here, f'{g) = 1/ I3{g). This is equivalent to the Gell-Mann — Low — Lee solution (plsl). 
To every solution of eq.(pJ|) for a function s, there corresponds some function S(?/, g) 
related to s by 

Y.[y/x,g{x,y]g)] 



s{x,y;g) 



(39) 



^iy,g) 

Let us give also the general solution of the same type for the system (^) for the k- 
couplings case. It can be written down in terms of k arbitrary functions Fi, reversible 
simultaneously with respect to last arguments, and defined from the system of k func- 
tional relations 



Fiiy,{g})=F, 



y 



X 



{gix,y; {g})} 



, {g} = gi,-- gk ; ij = i,-- k . (40) 



All solutions (|37| ) — (^) satisfy the usual normalization conditions. 
The transition to the massless limit in expressions (^) 
a trick analogous to the given above. Then, e.g.. 



(^) can be performed by 



s{x,g) = X 



^^{gix,g)} 
^9) 



Let us also formulate solution for the 2-coupling case gi = g, g2 = h in the massless 
limit in the form analogous to (Bl) 



fi{g,h) = fi{g,h) + lnx, 



1,2 . 



(41) 



From the solutions presented it follows that imposing group symmetry one reduces 
by unity the number of independent arguments. 



2.4 RGM algorithm 

2.4.1 Technology of RG Method 

The idea of the approximate solution marriage |^, ^ with group symmetry can be realised 
with help of group DEs. If we define group generators (3, 7 from some approximate 
solutions and then solve evolutional DEs, we obtain RG improved solutions that obey 
the group symmetry and correspond to the approximate solutions used as an input. 

Now we can formulate an algorithm of improving an approximate solution. The 
procedure is given by the following recipe which we illustrate by a massless one-coupling 
case (I) and (|): 

Assume some approximate solution ^appr is known. 

1. On the basis of eq.(22a) define the beta-function 

(42) 
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2. Integrate eq.(|), i.e., construct the function 

def 



(43) 



3. Resolve the eq.(p8D 

g^G{^^g) = r'{f{9) + ^^x} . (44) 

4. Then, the solution g^Q, precisely satisfies the RG symmetry, i.e., it is an exact 
solution of eq. and corresponds to ^g^ppr • 

For illustration, take as a ^appr ^^e simplest perturbative expression (pQ]) for the 
invariant coupling. Here, the /5-function is P{g) = —Pig'^ , and the integration yields 



9 dg _ 1 fl r 



The solution obtained 

g{x,g) = — — ^-j — , (45) 
1 + gpi In X 

one one hand, exactly satisfies the RG symmetry and, on the other, being expanded in 
powers of g, correlates with the input 



2.4.2 RGM usage in QFT 

As it has been explained above in Section the QFT perturbation expression of 
finite order does not obey the RG symmetry. On the other hand, in Section p.4.1| it 
was shown that the one- loop UV approximation for g used as an input in eq.(^) for 
the construction of a group generator P{g) yields expression (|45D that obeys the group 
symmetry and exactly satisfies FE (^). 

Now, using the geometric progression ( P5| ) as a hint, let us represent the 2-loop 
perturbative approximation for g in the form 

%t=g- g^Pi Inx + g^ [Pi In' X - /32 In a;] + 0{g^) , 

where /3i and ^2 mean the /? -function coefficients at the one-loop and two-loop level, 
respectively. If we substitute this expression into eq.(|^) we obtain 

g?hx,g)- g^^ [x/t- {t,g)) = g'Pl ln(x/t) In' t . 

Meanwhile, we can use 'g'^p^ as an input in Eq. (^2D . Now the step 1 yields 

P^'\g) = -Pig' - P2g' 

and then (step 2) 

^' J 72 + 673 l + hz' i3^- ^ ' 
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To make the last step, we have to start with the equation 



which is a transcendental one and has no simple explicit solution^. Due to this, one has 
to resolve this relation approximately. Take into account that the second, logarithmic, 
contribution to f'^'^\z) in (^61) is a small correction to the first one at hz <^ 1. Under 
this reservation we can substitute the one-loop RG expression ( ^5D instead of g^^ into 
this correction and obtain the explicit expression 

77(2) = ^ ■ / = In T (47) 

l + ^/5i/ + ^(/52//501n [1 + ^7^/] ' ■ ^ ^ 

This result (first obtained |^ in mid-50s) is interesting in several aspects. 

First, being expanded in g and gl powers, it produces an infinite series containing 
"leading", i.e. ~ g{gl)"', and "next-to-leading" ~ g'^{gl)^ UV logarithmic contributions. 
Second, it contains a nontrivial analytic dependence 

ln(l + (7/?iO ~ln(lng') 

which is absent the in perturbation input. Third, being compared with eg. (^51) , it 
demonstrates algorithm of subsequent improving of accuracy, i.e., of RGM regularity. 

Now we can resume the RGM properties. The RGM is a regular procedure of combin- 
ing dynamical information (taken from an approximate solution) with the RG symmetry. 
The essence of RGM is the following: 

1) The mathematical tool used in RGM is Lie differential equations. 

2) The key element of RGM is possibility of (approximate) determination of group 
generators from dynamics. 

3) The RGM works effectively in the case when a solution has a singular behaviour. It 
restores the structure of singularity compatible with RG symmetry. 



3 RG in QFT 

This section is devoted, mainly, to general topics of RG applications in the QFT short- 
distance asymptotic behavior. We discuss the specific features of UV analysis connected 
with use of perturbation theory, in particular, reliability of results. 



3.1 UV analysis in general 



•^It can be expressed in terms of a special, Lambert, W^-function : W{z) exp^^^^ = z; see, e.g., p9[ . 
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3.1.1 One-coupling case 

General analysis of the UV asymptotic behavior for the one-coupling QFT model can 
be performed rather simply on the basis of the solution 

of the massless RG equation (H) for an effective coupling with g = g{l,g) and /x, a 
reference point. As follows from it, the asymptotics at Inx ^ oo corresponds to the 
divergence at the upper limit of the l.h.s. integral. Depending on the feature of the 
/3-function the resultant UV behaviour of the invariant coupling g differs very much. 

Suppose that at very small g values the beta-function is positive. Then, three cases 
are possible : 

a) Consider first the situation with 

dz 

oo (52a) 



f3{z 

corresponding to the case when the beta-function has a zero at some finite point g^ . 

Here, the UV asymptotic value of effective coupling is finite g{oo,g) = g^: < oo , 
which relates to the finite renormalization of the coupling constant: Z = g{oo, g)/g. 

Using the terminology of DEs qualitative theory, one can say that dX g = g^ we have 
a UV fixed point. 

li aX g = (7* there is a first order zero [3{g) ~ h{g^: — g) , then eq.(52a) gives 

g{x, g)-g,^C exp^^^"^ = C{Q^/^l'^)-^ as ^ oo , (50) 

i.e., in the vicinity of a fixed point we have an asymptotic power legime. 

b) If (3{g) is (monotonically) increasing as (7 ^ oo but gentler than (7^, so that 

dz . , X 

W) ' °° 

then the effective coupling tends to infinity 

lim g{x,g) 00 , 

which corresponds to infinite coupling constant renormalization. Formally, this is equiv- 
alent to (7oo = 00. 

c) At 

d^ 



L = In a; 00 < 00 , (52c) 



that happens if 

lim P{g)/g'^ > const 
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the theory has an inner contradiction, as far as 

g{xoo,g) = 00 at Xoo < 00 

and the momentum region x > Xoo can not be described by the theory. We encounter 
here a ghost trouble, as explained below in this Section. 

Up to now we have assumed that the generator i3{g) is positive. In the opposite case 



d) i3{g) = —b{g) < one has to deal with the equation 

r9 dz 



biz) 



\nx {52d) 



and study possible divergence of the integral involved at the lower limit. 

If this occurs at some finite value g = g^o , the situation is quite analogous to the 
case b) . The only difference is that now the effective coupling tends to its limiting value 
goo from above. 

As the most important case, we consider the possibility when the singularity lies at 
the origin j3{0) = which happens in QCD. Then, g vanishes g{oo,g) = in the UV 
limit which corresponds to the asymptotic freedom phenomenon. E.g., if we assume here 
that i3{g) = —Pig"^ at g ^ , then 

9{x,9)^-T^ — X 00 . (51) 

Pi Inx 

3.1.2 Multi-coupling case 

For the quantum field model with several coupling constants one has to consider the 
system of coupled functional (^6[) or differential equations. The last ones can be analyzed 
by the well-known methods of the qualitative theory of differential equations. 

Take the case with two coupling constants g and h. The system of evolutional 
differential equations is 

^g = Pg{g,h), h = [3h{g,h)- f = df/di; i = \nx. (52) 

According to (^T]), the general solution to this system is of the form 

F{g, h) = F{g, h) + i , $(^, h) = <l>{g, h) + i . 

where F and $ — two arbitrary reversible functions. 

As far as argument I = lnx does not enter explicitly into the generators jSg and jSh, 
it can formally be excluded by dividing one of the equations (|5^) by the other : 

|=m^), F = ^. (53) 
dh fJh 

This equation can be analyzed on the two dimensional phase plane {g, h). 
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First explicit example of such phase portrait has been obtained in mid-fifties by 
I. Ginzburg [|1^ — see also Section 51.4 in the third edition of the monograph |jl|. 



The essential features are now singular points and singular solutions. Singular points 
correspond to /?j = (or = oo). They can be of different types: a stable fixed point that 
is known as attractor, an unstable fixed point and a saddle- type point. In the vicinity 
of the UV attractor one can have a power scaling behavior as in eq.(|SD|). Singular 
solution, separatrix, joins singular points and can also be stable or unstable. Generally, 
the unstable ones separate the parts of phase plane with different UV asymptotes that 
correspond to UV stable separatrices. 

3.2 Perturbative approach to the UV asymptote 

3.2.1 Structure of RG results 

Consider a general situation with the RG approach to the UV asymptotic behavior based 
on perturbation calculation input. In the one-coupling QFT case, group generators 
entering into DEs can be written as 

m = f3i9' + f329' + ... , ^{g) = Ag + ^29' + .... (54) 
Generally, expansion coefficients depend on the mass variable 

P{y,g) = J2My)9'''' , ^{y,g) = J2My)9' ■ (55) 
i>i I 

Note that if g is just the ^-matrix expansion parameter (that can be equal to the 
coupling constant or to its square) then usually the first term in expansion for f3 is 
quadratic and for ip — linear, as it is explicitly indicated above. 

Substituting (Q) into (^8]), and re-expanding the ratio 7^//?(7), we obtain after 
integration 

1/g -i/g-^ Hg/g) - Ug -g) + g') = A inx , (56) 

In s{x, g) = i^i/Pi) ln{g{x, g)/g} + C2{g - g) + 0(/) , (57) 
h = 13^1 13i - {I32/I3i)\ C2 = Vp2/l3i]{ip2/^i - /?2M). 

As follows, the solutions g and s depend on two arguments g and g\\ix . By expanding 
them in powers of g we get 

g{.x,g) = gfi{g\nx)+g'^ f2{g\Tix) + . . . , \Yis{x,g) = ipi{g\Y).x)+gip2{g\Tix) + . . . (58) 

where and ipi have a simple form. For example, fi{z) = (1 — [3izy^ , '{'li.z) ~ /2(^) ~ 
ln/i(^) . 

Comparing expressions obtained with usual perturbative expansions 
gpt{x, g)=g + g^Pi Inx + g^[Pl Inx + /?2 Inx] + 0(/), 
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Spt{x,g) = l + ^V^ilnx + /[(/5iV^i/2)ln2x + ^2lnx] +0(^^) , 

used as an input to obtain our starting generators, one can see the qualitative effect 
of the RGM using. In the case considered, it changes the region of applicabihty of the 
perturbation method limited by the condition In x <^ 1 to a more larger region defined 
by two relations 

g{x,g)<t:l, (59) 

the second of which is defining. 



3.2.2 The ghost-pole trouble 

Turn now to the one-loop RG approximation for the effective coupling g, considered in 
the UV, i.e., massless limit. According to (^) and (p3D, it has the form 



9{i){x,g) ^ 



1 — Pig In X 

Let the numerical coefficient /3i be positive. Such is the case in QED where Pi = l/3n 
and g stands for the expansion parameter a = e^. This expression obviously has a pole 
singularity at 

X = x^ = exp{l/Pig) = exp(37r/a) . 

As far as the QED effective coupling is proportional to the (transverse part of) photon 
propagator, this pole, generally, describes some bound state of a system with the photon 
quantum numbers. However, a pole related to a physical bound state must have positive 



residue while the l.h.s. of eq.(23) has a negative one. 

This means that it corresponds not to a physical but rather to some unphysical, 
so-called ghost, state. The presence of a ghost singularity can be treated as a signal of 



inconsistency of a theory. Such claims have been made pO[ in the mid 50s when the 



ghost-pole trouble was first discovered [|2T| just before the birth of the RGM. 

The RG method proved to be very effective for a general discussion of the ghost-pole 
issue. The first question that must be answered here is the stability of indication of the 
ghost-pole existence with respect to the multi-loop corrections. 

Note that in a perturbation calculation, the /5-function depends on the adopted 
renormalization procedure; at the massless case, starting with the 3-loop level the coef- 
ficients of the perturbation series ( ^^ depend on the renormalization scheme (RS) used. 
In QED, the 3-loop P function in MOM (i.e., momentum subtraction) scheme is 

''--i^-£;-^(^c(3)-^). («o) 

The numerical value of the last parenthesis is about 0.4. Neglecting it for the moment, 
we start our discussion with the two-loop approximation for the P function. According 
to eq.(0), the 2- loop iterative RG solution is 

a 

"(2) = 1 a, , 3a 



3tt 
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This solution has an error of an order of a'^l and is interesting from several points 
of view. As it has been mentioned before, its a expansion besides leading logs contains 
an infinite number of next-to- leading terms a'^{al)"^, the first of which has been used 
as an input for construction of the j3 function. Second, in the vicinity of the ghost pole 
of the one-loop RG solution at li = 3n/a, the two-loop a(2) solution differs from 
considerably. Hence, an infinite sum of the next-to-leading logarithmic contributions in 
the region a/ ~ 1 becomes important. 

It is not trivial because for an each order of the perturbation input the next-to-leading 
term is negligible comparing with the leading one of the previous order (the ratio being 
of an order of aP2/Pi = Sa/An ~ 2.10"^). It can be seen with the help of ( ^6|) that the 
allowing for the last, 3-loop, term in (|60|) also becomes essential for the a ~ 1 case. 

This means that the problem of existence for the ghost pole in QED cannot be 
solved by taking into the account of next-to-leading and so on logs. Moreover, one can 
argue p2| on general RG ground that it is impossible to make any qualitative statement 
about the UV asymptote for the P{g) > g case basing on RG-improved perturbation 
calculations. Our next example illustrates this thesis. 



3.2.3 Scalar quartic model 

For the nonlinear scalar field with the quartic (self) interaction Lagrangian 



47r2 



int 



important progress has been achieved in 80s in the higher perturbation orders calcula- 



tion. The P function was calculated in the MS-scheme up to the 5- loop level |23 



MS 



-.9 



17 
"6 



■g^ + 16.275^' - 135.8^^ + 1437^^ 



We see from this expression that, in conrast with the QED case, due to its alternate- 
sign structure, there is no stability here even on a qualitative level. The odd-order 
approximations have a ghost-pole type behavior, whereas the even ones yield the fixed 
point (finite charge renormalization) case. Note also that, as can be shown [^], the 
upper boundary of the 10% confidence region corresponds to g values close to 0.1 . 

To comprehend the "loop dependence" of this boundary it is useful to represent 
expression f^j^^^ in a slightly different form 



r'i9) = ^9' 



1 - 



9 



0.529 



+ 



9 



0.303 



— )' 

0.222/ 



9 



0.180 



4n 



(61) 



It is clear now, that a boundary of the confidence region diminishes with rising the 
order of a loop approximation. The expression ( pl| ) looks like a beginning of a power 
asymptotic series of the Poincare type. Indeed, if we represent the f3 -function in a series 
expansion form (K3) then, as it is can be shown, the coefficients /?„ at n ^ 1 behave 
like ~ n\. 
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The method of determining asymptotic estimates of the perturbation expansion co- 
efficients of the Green functions uses a representation in the form of a functional (i.e. 
path) integral. This integral written down for the mentioned expansion coefficient can 
be calculated by the steepest descent method in the function space. To the saddle point, 
there corresponds an "instanton"-type Euclidean classical solution with a finite action. 

In this manner, an asymptotic expression was obtained E^ for the coefficients of the 
P function expansion. It has the form 

/3„ ^ n'l^ (l + 0[n^)\ (n ^ oo) . (62) 

The factorial growth of coefficients indicates that this is a power asymptotic series with 
zero radius of convergence that cannot be summed in the usual manner. We can obtain 
the information about the singularity structure at the origin {^g = 0) by using some 
special procedures. One of them is the Borel summation method. 

Here, we give short exposition of the results on the attempt of the summation of the 
series of (|6T|) type made in |Q (see also the review p6| ). 



Authors of the Ref. used as an input the /5-function 4-loop expression in the 
symmetric MOM-scheme 

pr^'ia) = - ^^g' + 19.3^^ - 146/ . (63) 

This alternate-sign asymptotic series can be summed by the Borel method. The idea 
is to represent the sum in a form of a Laplace transform integral. It is not difficult to 
see that the transition to the Laplace image just "kills" the factorial factor n\. For the 
modified Borel transformation 



-exp{-x/g)[x—j B{x) (64) 



dx , , , ( d 

g 

perturbation series can be written down as 



B{x) = Y^^x^ 



It has a nonzero circle of convergence and can be summed within the circle. However, as 



the integration domain in (^4|) goes outside the convergency region, we must make an an- 
alytic continuation for the function B[x). It can be done by a conformal transformation 
of the x-plane into the w-plane to map the domain of integration [0, oo] into the interior 
of the unit disk and the cut [—00, —1] into the boundary of the disk. One can choose 
this transformation in such a way that it correctly reproduces the singularity on the cut. 
The result of conformal transformation x ^ w{x) = (vT+~x — 1) ■ (vT+~r + l)~^ "looks 
quite well" : 

3x2 

B{x) = —(1 - 0.32w; - 0.127^2 + 0.084u7^) . 
128 
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Then, by transformation reverse to ( |5^ one can reconstruct beta function (3{g) which 
is nonanalytic in the g variable with essential singularity at g = 0. Graphs of the 
function (3{g) obtained by the Borel summation with allowance for the 1-, 2-, 3- and 
4-loop approximations look very similar one to other. 

They all lie now in a narrow parabolic ray slightly below the original one-loop 
parabola and within the limit of 10% accuracy enable to advance into the region g ~ 50. 
This means that the summation procedure adopted enlarges the confidence interval in 
several hundred times! Besides this it gives the qualitative stability of results. All they 
are now in favour of a ghost-type UV asymptote. 

Nevertheless, these results can be considered only as a support but not the proof 
of the 0^-model inconsistency. The weak point here is that starting with (|6^) we have 
assumed the definite analyticity properties of P{g) in the whole complex (yf-plane. 



3.3 Mass— dependent analytic solution 

A general method of an approximate solution to the massive (i.e., mass-dependent) RG 
equations was developed in Ref . . Analytic expressions of a high level of accuracy for 
an effective coupling and a one-argument function were obtained up to 3- and 4-loop 
order |^ . 



For example, the two-loop massive RG-solution for the invariant coupling 



2n(2) 



1 + asAiiQ^, m?) + a.AsMi) In (1 + 
at small as values corresponds to perturbation expansion 



(65) 



«.(Q') 



2^(2) 
pert 



a. 



+ ... . 



At the same time, it smoothly interpolates between two massless limits (with Ai ~ 
jSeliaQ^ + q) at <^ and ^ described by an equation analogous to (^T]). 
In the latter case it can be represented in the form usual for the QCD practice: 

a;\QVA')i^ ^ jln § + &i In (in fj) } ; 6i = | • 

Solution ( |65D demonstrates, in particular, that the threshold crossing generally chan- 



ges the subtraction scheme |p9 
The investigation [[ 



28 1 was prompted by the problem of taking explicitly into 



account of heavy quark masses in QCD. However, the results obtained are important 
from a more general point of view for a discussion of the scheme dependence problem 
in QFT. The method used could also be of interest for RG applications in other fields 
within the situation with disturbed homogeneity, such as, e.g., intermediate asymptotics 
in hydrodynamics, finite-size scaling in critical phenomena and the excluded volume 
problem in polymer theory. 
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In paper this method was apphed to the evolution of effective gauge couphngs 
in Standard Model (SM). Here, a new analytic solution of a coupled system of three 
mass-dependent two-loop RG equations for three SM gauge couplings was obtained. 

For this goal, one has to start with a perturbative input for the SM couplings 



(Q^m^) ~ai-a2A(Q,m,/i) + afA2(Q)-a,'^a,Ai,(Q) ; z = l,2,3. (66) 



where Ai and Aij are one- and two-loop mass- and /x-dependent contributions of ap- 
propriate Feynman diagrams. In the framework of a massive renorm-group formalism 
1^, H, [ly] the corresponding Lie equations look like 



d,(Q^m^)~-a/(Q) 



A{Q) + Y.a,{Q)A,{Q) 



(67) 



with A = dA/di; i = InQV/x^ . Note that in the UV limit Ai{Q) = ; (Q) = (3, 
we arrive at the system 



that is commonly used - see, e.g., Refs. ||3T| - for the discussion of data extrapolation 
across the gauge desert and possibility of Grand Unification. 
The latter system can be solved iteratively in the form 



(68) 



Here, we present a generalization of this solution for the massive case, that is conve- 
nient for taking into account of threshold effects and discussing, in particular, the issue 
of the Grand Unification consistency check. 



Using the method of the paper one can obtain explicit iterative solution to the 
system (^7[). Here, as in the massless case, one first solves the one-loop approximation 

to (§3) to get g 

ai'\Q',m^) = [l/a, + A{Q^,m^)]-' . 

Inserting then this explicit expression into the second factor in the r.h.s. of (^) and 
performing an approximate integration of some integral — for detail see paper 
we arrive at the expression 



aj(Q^m2) ai 
quite analogous to (|65|). 



^ +^.(Q) + e4#Mi + «A(Q)]- 



r MQ) 



(69) 



This exact solution of an one- loop massive RG equation was first obtained in Ref. |3J 
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The remarkable feature of this solution is that it depends explicitly only on mass- 
dependent perturbation coefficients Ai{Q), Aij{Q) and, being expanded in powers of 
coupling constants, exactly corresponds to the perturbative input (|66|) . On the other 
hand, in the massless limit it goes to solution (p^). 

The accuracy of the last approximate expression can be estimated by the method 
used in paper [^] . Generally, it corresponds to the accuracy Q of three-loop expression 
( ~ In ) for the effective coupling in the one-coupling case that is quite sufficient for 
current pactice. 



3.4 Some important results 

In the early 70s, S. Weinberg proposed the notion of a running mass of a fermion. 



If considered from the viewpoint of paper [|T3[, this idea can be formulated as follows: 

any parameter of the Lagrangian can be treated as a (generalized) coupling constant, 
and its effective counterpart should be included into the renorm-group formalism. 

New possibilities for applying the RG method were discovered when the technique of 
operator expansion at short distances (on the light cone) appeared The plausibility 



of this approach stems from the fact that the RG transformation, regarded as a Dyson 
transformation of the renormalized vertex function, involves the simultaneous scaling 
of all its invariant arguments - normally, the squares of the momenta. Meanwhile, for 
the physical amplitude, some of them are fixed on a mass shell. The expansion on the 
light cone, so to say, "separates the arguments" , as a result of which it becomes possible 
to study the physical UV asymptotic behaviour by means of the expansion coefficients 
(when some momenta being fixed on mass shell). As an important example, we can 



mention the evolution equations for moments of QCD structure functions |^ . 

The revealing of an asymptotic freedom phenomenon can be considered as the most 
important result obtained in particle physics by the RG technique. 



Historically, this discovery was made in the framework of the SU(3) non-Abelian 



Yang-Mills model in the early 70s. Since that time this model for the eight-component 
4- vector field B'^^{x) was adopted as a basic ingredient for the QFT description of matter 
on the parton level. 

The key point is that self-interaction of this non-abelian gluonic quantum field due to 
dominance of its unphysical components gives negative contribution to the beta function 
perturbation expansion. For the two- loop (scheme-independent) case 

with positive /3i^2 for a number of quark fiavours small enough. 
Correspondingly, the one-loop renorm-group expression 

-(1)/ N '^■'^ 



1 + asPi Inx ' 



^See eqs.(13) and (16) in Ref.KJ 
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for the QCD effective coupling exhibits a remarkable UV asymptotic behaviour thanks 
to Pi being positive. This expression implies, in contrast to eg. (j23|) , that the effective 
QCD coupling decreases as x ~ increases and tends to zero in the UV limit. This 
feature discovered in the early 70s, precisely corresponded to the parton physical picture 
of the hadronic structure. 

One more interesting application of the RG method in the multicoupling case, as- 
cending to 50s ||12[ , refers to special solutions, the so-called separatrixes in a phase space 
of several invariant couplings. These solutions relate effective couplings and represent 
scale-invariant trajectories, like, e.g., gi = gi{gi) in the phase space which are straight 
lines in the one-loop case. 

Some of them that are "attractive" (or stable) in the UV limit, are related to sym- 
metries that reveal themselves in the high-energy domain. It was conjectured that these 
trajectories may be related to hidden symmetries of a Lagrangian and even could serve 
as a tool to find them. On this basis the method was developed ||39| for finding out these 



symmetries. It was shown that in the phase space of invariant couplings the internal 
symmetry corresponds to a singular solution that remains a straight-line when higher or- 
der corrections are taken into account. Such solutions corresponding to supersymmetry 
were derived for some combinations of gauge, Yukawa and quartic interactions. 
Generally, these singular solutions obey the relations 

dgi _ dgi_dgi_ ^ _ 
dl dgi dl ' 

which are known since Zimmermann's paper as the reduction equations. In the 80s 
they were used |^ (see also review paper |^2| and references therein) in the UV analyzis 



of asymptotically free models. Just for these cases the one-loop reduction relations are 
adequate to physics. 

Quite recently some other application of this technique was obtained in supersym- 
metric generalizations of Grand Unification scenario in the Standard Model. It was 
shown p3| , Q that it is possible to achieve complete UV finiteness of a theory if Yukawa 



couplings are related to the gauge ones in a way corresponding to these special solutions, 
that is, to reduction relations. 



The mass- dependent technique described in Section O] was successfully used for 



the development of the Dhar-Gupta approach |]45| , |46[] that led to finite perturbative 



predictions for a physical quantity which is free of renormalization scheme ambiguities. 



In paper [47|, this approach was reformulated for the mass-dependent case with several 



coupling constants. 

One more recent QFT development relevant to the renorm- group is the "Analytic 
approach" to perturbative QCD (pQCD). It is based upon the procedure of Invariant 
Analyticization ^ ascending to the end of 50s. 

The approach consists in the combining of two ideas: the RG summation of UV logs 
with analyticity in the variable, imposed by spectral representation of the Kallen- 
Lehmann type which implements general properties of the local QFT including the 
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Bogoliubov condition of microscopic causality. This combination was first devised 
to get rid of the ghost pole in QED about forty years ago. 

Here, the pQCD invariant coupling asiQ"^) is transformed into an "analytic cou- 
phng" Q;an(<5^/A^) = A{x) which, by construction, is free of ghost singularities due to 
incorporating some nonperturbative structures. 

This analytic coupling A{x) has no unphysical singularities in the complex Q^-plane; 
its conventional perturbative expansion precisely coincides with the usual perturbation 
one for asiQ"^) ] it has no extra parameters; it obeys a universal IR limiting value ^(0) = 
4:n/Po independent of the scale parameter A; it turns out to be remarkably stable 
in the IR domain with respect to higher-loop corrections and, in turn, to the scheme 
dependence. 

Meanwhile, the "analyticized" perturbation expansion |5l| for an observable F, in 
contrast to the usual case, may contain specific functions An{x), instead of powers 
{A{x))"' . In other words, the pertubation series for F{x), due to analyticity imperative, 
can change its form in the IR region turning into an asymptotic expansion a la 
Erdelyi over a nonpower set {An{x)} . 



4 RG expansion 

In 70s and 80s RGM was applied to (besides QFT) critical phenomena: poljTiiers, tur- 
bulence, non-coherent radiation transfer, dynamical chaos, and so on. Simpler and less 
sohpisticated motivation in critical phenomena (than in QFT) makes this "explosion" 
of RG applications possible. 



4.1 Critical phenomena 

4.1.1 Spin lattices 

The so called renormalization group in critical phenomena is based on the Kadanoff- 



Wilson procedure |5^, |5^ referred to as "decimation" or "blocking" . Initially, it emerged 
from the problem of spin lattice. Imagine a regular (two- or three- dimentional) lattice 
consisting of A^"', d = 2,3 cites with an 'elementary step' a between them. Suppose, that 
at every site a spin vector a is sitting. The Hamiltonian describing the spin interaction 
of nearest neighbours 

H = kJ2a,-a,±, (70) 

i 

contains k, the coupling constant. The statistical sum is obtained from the partition 
function, 5* =< exp{—H/9) >aver. ■ 

To realize the blocking or decimation, one has to perform the "spin averaging" over 
block consisting of n*^ elementary sites. This is a very essential step as far as it dimin- 
ishes the degree of freedom number (from N'^ to (N/nY). It destroys the small-range 
properties of the system under consideration, in the averaging course some information 
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being lost. However, the long-range physics (like correlation length essential for phase 
transition) is not affected by it, and we gain simplification of our problem. 

After this procedure, new effective spins S arise in sites of a new effective lattice 
with a step na. We obtain also a new effective Hamiltonian, with new effective coupling 
Kn that has to be defined in the averaging process as a function of k and n 



where AH contains quartic and higher terms; AH = Yl^Ti ■ "S + 

For the IR (long-distance) properties, AH is unessential. Hence, we can conclude 
that the spin averaging leads to an approximate transformation. 



a ■ a 



(71) 



or, taking into account the "elementary step" change, to 



KWn-. {a 



na, 



The latter is the Kadanoff- Wilson transformation. 

In general, the "new" coupling constant Kn is a function of the "old" one and of 
the decimation index n. It is convenient to write it down in the form Kn = K{l/n, K). 
Then, the KW transformation can be formulated as follows: 



na, 



k^Kn 



K 



KW{n) 

These transformations obey the group composition law 



n 



(72) 



KW{n) ■ KW{m) = KW{nm) 



if 



K{x,k) 



,x 



K{-,K{t,k)) 



X 



t 



nm 



(73) 



This is just the RG symmetry. 

We observe the following points: 

• The RG symmetry is approximate (due to neglecting by AH). 

• The transformations KW{n) are discrete. 

• There exist no reverse transformation to KW{n) . 

Hence, the 'Kadanoff- Wilson renormalization group' is an approximate and discrete 
semi-group. For a long distance (IR limit) physics, however, AH is irrelevant, A{l/n) 
is close to continuum and it is possible to use differential Lie equations. 

In application of these transformations to critical phenomena the notion of a fixed 



point is important. As it was explained in Section |3.1| , it is usually associated with 
power- type asymptotic behavior. Note here that, contrary to the QFT case considered 
in Section RTT], in phase transition physics we deal with the IR stable point. 
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4.1.2 Polymer theory 

In the polymer physics one considers statistical properties of polymer macromolecules 
which can be imagined as a very long chain of identical elements. The number of elements 
N could be as big as 10^, the macromolecular size reaching several hundred Angstroms. 

Such a big molecular chain forms a specific pattern resembling the pattern of a 
random walk. The central problem of the polymer theory is very close to that of a 
random walk and can be formulated as follows. 

For a very long chain of N "steps" (the size of each step = a) one has to find 
the "chain size" i?Ar as the distance between the "start" and the "finish" points, the 
distribution function of angles (pi between neighboring elements being given. 

The function /(0) is defined by the forces between adjacent elements depending on 
some external factors like temperature T. The essential feature of a polymer chain is 
the impossibility of a self-intersection. This is known as an excluded volume effect in 
the random walk problem. In reality, polymer molecules are swimming in a solvent and 
form globulars. 

For large values the molecular size Rn follows the power Fleury law Rn ~ A^"^ 
with i/, the Fleury index. When A^ is given, i? at is a functional of /(</>) which depends on 
external conditions (e.g., temperature T, properties of solvent, etc. ). If T increases, Rn 
increases and at some moment globulars touch one another. This is the polymerization 
process very similar to a phase transition phenomenon. 

The Kadanov-Wilson RG (KWRG) blocking ideology has been used in polymer 
physics by De Gennes The key idea is a grouping of n chain subsequent elements 



into a new "elementary block". This grouping operation is very close to Kadanoff's 
blocking. It leads to the transformation 

{1 ^ n ] a An} 

which is analogous to one for spin lattice decimation. This transformation must be 
specified by a direct calculation which gives the explicit form of An = d{n,a). Here we 
have a discrete semi-group. Then, by using the KWRG technique, one finds the fixed 
point, obtains the Fleury power law and can calculate its index u. 

Generally, the excluded volume effect yields some complications. However, inside the 
QFT RG framework it can be treated rather simply |^ by introducing an additional 
argument similar to finite length L in transfer problem and particle mass m in QTF. 

Besides polymers, the KWRG approach has been used in some fields of physics, like 
percolation, noncoherent radiation transfer ||56| , dynamical chaos ||5^ and some others. 



Meanwhile, the original QFT-RG approach proliferated into the theory of turbulence. 
4.1.3 Turbulence 

To formulate the turbulence problem on the "RG language" one has to perform the 
following steps pUI, pO, BT|: 



1. Introduce the generating functional for correlation functions. 
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2. Write the path integral representation for this functional. 

3. By changing the function integration variable find the equivalence of the classical 
statistical system to some quantum field theory model. 

4. Construct the system of Schwinger-Dyson equations for this equivalent QFT. 

5. Perform the finite renormalization procedure. 

6. Derive the RG equations. 



4.2 Paths of RG expansion 

RG is expanded in diverse fields of physics in two different ways: 

• by direct analogy with the Kadanov- Wilson construction (averaging over some set 
of degrees of freedom) in polymers, non-coherent transfer and percolation, i.e., 
constructing a set of models for a given physical problem. 

• search for an exact RG symmetry by proof of the equivalence with a QFT model: 
e.g., in turbulence (Refs. §T|), plasma turbulence and some others. 



To the question Are there different renormalization groups? the answer is positive: 

1. In QFT and some simple macroscopic examples (like, one-dimentional transfer 
problem) , RG symmetry is an exact symmetry of the solution formulated in its 
natural variables. 

2. In turbulence, continuous spin-field models and some others, it is a symmetry of 
an equivalent QFT model. 

3. In polymers, percolation, etc. , (with KW blocking), the RG transformation is 
a transformation between different auxiliary models (specially constructed for this 
purpose) of a given system. 

As we have shown, there is no essential difference in the mathematical formalism. 
There exists, however, a profound difference in physics: 

— In cases 1 and 2 (as well as in some macroscopic examples), the RG is an exact 
symmetry of a solution. 

— In the Kadanov- Wilson-type problem (spin lattice, polymers, etc. ), one has to 
construct a set Ai of models Mj. The KWRG transformation 

R{n)Mi = Mni , with integer n (74) 

is acting inside a set of models. 
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4.3 Two faces of RG in QFT 

As it was explained in Section 1.4, the vacuum, i.e., the interparticle space, contain 
vacuum fluctuations. Due to them, the charge of a particle is screened. In accordance 
with Dirac eq.(^), in momentum space the dependence of an electron charge can be 
presented 

e(g2)=e|l + ^ln(QV,2) + ...} ; = e^(l/r,2) = 1/137. (75) 

in terms of the classical electron charge and of electron Compton lenght. 

The first idea of an additional symmetry in this problem was born by Stiieckelberg 
and Peterman [^] . In their pioneering investigation the very existence of group transfor- 
mation was discovered within the renormalization procedure the result of which contains 
finite arbitrariness. Just this degree of freedom in finite renormalized expressions was 
used by Bogoliubov and Shirkov in Refs.[0] — []TU|. Roughly speaking, this corresponds 
to the change re — > l//i. 

The basic idea was that, instead of l/rg, one can use some other reference point [i. 
This is equivalent to introducing of a new degree of freedom associated with the reference 
point scale. Instead of (^) we have 



e{QVf.') = e,\l + ^ln^ + ...\ . (76) 



Here, the effective charge is considered after the subtracting of infinities and is given 
by a "finite representation" (|76[). The RG symmetry is formulated in terms of the 
scale and /i represents the reference point. 

Another approach was used by Gell-Mann and Low 0. Their paper was devoted to 
the short distance behaviour in a nonlocal QED with a cutoff A, and the "A degree of 
freedom" was used to analyze the UV behaviour. Instead of renormalization, there is a 
regularization and the charge is given by the "singular representation" 

e(A)=efl + ;^lnAV2 + ... ] (77) 
V 67r / 

which is singular in the limit A ^ oo. 

We can draw a transparent picture (as was commented later by Wilson in his Nobel 
lecture) of the last approach. Imagine an electron of a finite size, smeared over a small 
volume with the radius Ri = h/cAi , ln(A^/mg) ^ 1. The electric charge Ci of such a 
non-local electron is considered as depending on the cut-off momentum Aj so that this 
dependence accumulates the vacuum polarization effects which, in reality, take place at 
distances from the point electron smaller than R^. We deal with a set of models of the 
non-local electron correponding to different values of the cut-off Aj. Here, depends on 
Ri and the vacuum polarization effects in the excluded volume R^ should be subtracted. 
In this language, the RG transformation is the transition from one value of the smearing 
radius to another Ri Rj, simultaneously with a corresponding change of the effective 
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electron charge Cj — >■ Cj. In other words, the RG symmetry here is that related to 
operations in the space of models of non-local QED constructed in such a way that at 
large distances every model is equivalent to the real local one. 



5 RG symmetry in mathematical physics 
5.1 Functional self-similarity 

The RG transformations discussed above have close connection with the concept of a 
self-similarity(SS). The SS transformations for problems formulated by nonlinear partial 
DEs are well known, since the last century, mainly in dynamics of liquids and gases. 
They are one parameter A transformations defined as simultaneous power scaling of 
independent variables z = {x,t, . . .} , solutions fk{z) and other functions Vi{z) 

S,: {x-.x\,t-^t\\ h{z) ^ fl{z') = X'^'^Mz') , V,iz) X-'^V.iz') } 

entering into the equations. 

To emphasize their power structure, we use a term power self-similarity = PSS. 



According to Zel'dovich and Barenblatt, [^, |6^ the PSS can be classified as: 

a/ PSS of the 1st kind 
with all indices a, ...y?, z/, ... being (half) integers (Integer PSS) that are usually found 
from the theory of dimensions; 

b/ PSS of the 2nd kind 
with irrational indices (Fractal PSS) which should be defined from dynamics. 

To relate RG with PSS, let us turn to the solution of the renorm-group FE 

9{xt,g) = g{x,g{t,g)) . 

Its general solution is known; it depends on an arbitrary function of one argument - 
see eg . (pSf) . However, at the moment we are interested in a special solution linear in 
the second argument: g{x,g) = gf{x). The function /(x) should satisfy the equation 
f{xt) = f{x)f(t) with the solution /(x) = x'^ . Hence, g{x,t) = gx^ . This means that in 
our special case, linear in g, the RG transformation (^ is reduced to PSS transformation, 

R^^{x^ xt-\ g ^ gf} = St . (78) 

Generally, in RG, instead of a power law, we have arbitrary functional dependence. 
Thus, one can consider transformations (P), ([T5| ) and (17) as functional generalizations 



of usual (i.e., power) self-similarity transformations. Hence, it is natural to refer to them 
as to the transformations of functional scaling or functional (self )similarity (FS) rather 
than to RG-transformations. In short, 

RG = FS , 

with FS standing for Functional Similarity. 



34 



We can now answer the question concerning the physical meaning of the symmetry 
underlying FS and the Bogoliubov's renorm-group. As we have mentioned, it is not 
a symmetry of the physical system or the equations of the problem at hand, but a 
symmetry of a solution considered as a function of the relevant physical variables and 
suitable boundary conditions. A symmetry like that can be related, in particular, to 
the invariance of a physical quantity described by this solution with respect to the way 
in which the boundary conditions are imposed. The changing of this way constitutes a 
group operation in the sense that the group composition law is related to the transitivity 
property of such changes. 

Homogeneity is an important feature of the physical systems under consideration. 
However, homogeneity can be violated in a discrete manner. Imagine that such a discrete 
inhomogeneity is connected with a certain value of x , say, x = y. In this case the RG 
transformation with the canonical parameter t will have the form (0) with the group 
composition law ([151) . 

The symmetry connected with FS is a very simple and frequently encountered prop- 
erty of physical phenomena. It can easily be "discovered" in numerous problems of 
theoretical physics like classical mechanics, transfer theory, classical hydrodynamics, 
and so on ^ 



5.2 Recent application to boundary value problem 

Recently, some interesting attempts have been made to use the RG concept in classical 
mathematical physics, in particular, to study strong nonlinear regimes and to investigate 
asymptotic behavior of physical systems described by nonlinear PDEs. 

About a decade ago, the RG ideas were applied by late Veniamin Pustovalov with 
co-authors to analyze a problem of generating higher harmonics in plasma. This 
problem, after some simplification, was reduced to a couple of partial DEs with the 
boundary parameter - "solution characteristic" - explicitly included. It was proved 
that corresponding solutions admitted an exact symmetry group that takes into account 
transformations of this boundary parameter, which is related to the amplitude of the 
magnetic field at a critical density point. The solution symmetry obtained was then 
used to evaluate the efficiency of harmonics generation in cold and hot plasma. The 
advantageous use of the RG-approach in solving the above particular problem gave 
promise that it may work in other cases and this was illustrated in |^ by a series of 
examples for various boundary value problems. 

Moreover, in Refs. [^, Q the possibility of devising a regular method for finding a 
special class of symmetries of solution to the boundary value problem (EVP) in math- 
ematical physics, namely, RG-type symmetries, was discussed. The latter are defined 
as solution symmetries with respect to transformations involving parameters that enter 
through the equations as well as through the boundary conditions in addition to (or 
even rather than) the natural variables of the equations. 

As it is well known, the aim of the modern group analysis [^, |7^, which goes back 



to works by S. Lie ||71|, is to find symmetries of DEs. This approach does not include a 
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similar problem of studying the symmetries of solutions of these equations. Outside the 
main direction of both the classical and modern analysis, there remains as well a study 
of solution symmetries with respect to transformations involving not only the variables 
present in the equations, but also parameters entering into the solutions from boundary 
conditions. 

From the afore-said it is clear that the symmetries which attracted attention in 
the 50s in connection with the discovery of the RG in QFT were those involving the 
parameters of the system in the group transformations. It is natural to refer to these 
symmetries related to FS (or RG-type) symmetries. 

It should be noted that the procedure of revealing the FS symmetry (FSS), or some 
group feature, similar to the FS regularity, in any partial case (QFT, spin lattice, poly- 
mers, turbulence and so on) up to now is not a regular one. In practice, it needs some 
imagination and atypical manipulation "invented" for every particular case — see the 
discussion in . By this reason, the possibility to find a regular approach to construct- 



ing FSS is of principal interest. 

Recently, a possible scheme of this kind was presented as applied to a mathematical 
model that is described by a BVP. The leading idea [^, |73[ in this case is based on 
the fact that solution symmetry for this system can be found in a regular manner by 
using the well-developed methods of modern group analysis. 

The scheme that describes devising of FSS and its application is then formulated 
jT^ , [76[| as follows. Firstly, a specific RG-manifold should be constructed. Secondly, some 
auxiliary symmetry, i.e., the most general symmetry group admitted by this manifold is 
to be found. Thirdly, this symmetry should be restricted on a particular solution to get 
the FSS. Fourthly, the FSS allows one to improve an approximate solution or, in some 
cases, to get an exact solution. 

Depending on both a mathematical model and boundary conditions, the first step of 
this procedure can be realized in different ways. In some cases, the desired FS-manifold 
is obtained by including parameters, entering into a solution via an equation(s) and a 
boundary condition, in the list of independent variables. The extension of the space of 
variables involved in group transformations, e.g., by taking into account the dependence 
of coordinates of the renorm-group operator upon differential and/or non-local variables 
(which leads to the Lie — Backlund and non-local transformation groups [^) can also 
be used for constructing the FS-manifold. The use of the Ambartsumian invariant 
embedding method and of differential constraints sometimes allows reformulations 
of a boundary condition in a form of additional DE(s) and enables one to construct 
the FS-manifold as a combination of original and embedding equations (or differential 
constraints) which are compatible with these equations. At last, of particular interest is 
the perturbation method of constructing the FS-manifold which is based on the presence 
of a small parameter. 

The second step, the calculating of a most general group Q admitted by the FS- 
manifold, is a standard procedure in the group analysis and has been described in detail 
in many texts and monographs - see, for example, |75|. 
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The symmetry group Q thus constructed cannot as yet be referred to as a renorm- 
group. In order to obtain this, the next, third step should be done which consists in 
restricting ^ on a solution of a boundary value problem. This procedure utilizes the in- 
variance condition and mathematically appears as a "combining" of different coordinates 
of group generators admitted by the FS-manifold. 

The final step, i.e., constructing analytic expression for the solution of the boundary 
value problem on the basis of the FS, usually presents no specific problems. 

A review of the results, which were obtained on the basis of the formulated scheme, 
can be found, for example, in [^, |7^, |80| . 

We mention briefly, the FS analysis result for a particular problem of nonlinear 
optics, the problem of the laser beam self-focusing in a nonlinear medium. Here, one 
have a EVP for a coupled system of two nonlinear PDEs with the boundary condition 
given in a form of two one-argument functions. With help of RG=FS approach one new 
exact analytic and one new approximate analytic solution (for the practically important 
Gaussian initial transverse profile) has been found 

The important qualitative features of this example are: 

- the two-dimension structure has been analysed, that has a singularity in the be- 
havior of derivatives with respect to transverse coordinate, 

- the algebraic structure of the FSS operators is different from that of "usual RG of 
the QFT type". Here, we meet with a set of several infinitesimal renormgroup oper- 
ators, each of those can be used to reconstruct the analytic (exact or approximate) 
solution of EVP starting from the perturbative theory result. Moreover, renorm- 
group operators for exact solutions of EVP obtained appear as Lie-Eacklund {not 
Lie) infinitesimal operators. 

Up to now the outlined regular method is feasible for systems that can be described 
by DEs and is based on the formalism of modern group analysis. However, it seems also 
possible to extend this approach to boundary value problems that are not described just 
by differential equations. A chance of such an extension is based on recent advances 
in group analysis of systems of integro-differential equations which allow transfor- 
mations of both dynamical variables and functionals of a solution to be formulated 



83[] . More intriguing is the issue of a possibility of constructing a regular approach for 
more complicated systems, in particular to those having an infinite number of degrees 
of freedom. The formers can be represented in a compact form by functional (or path) 
integrals. 
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